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Abstract 

For a graph G and a positive integer k, the graphical Stirling number 
S(G, k) is the number of partitions of the vertex set of G into k non-empty 
independent sets. Equivalently it is the number of proper colorings of G that 
use exactly k colors, with two colorings identified if they differ only on the 
names of the colors. If G is the empty graph on n vertices then S(G, k) reduces 
to S(n, k), the familiar Stirling number of the second kind. 

In this note we first consider Stirling numbers of forests. We show that 
if (Fn^) n >o is an y sequence of forests with having n vertices and 

c(n) = 0(^/72/ log n) components, and if Xn^ is a random variable that 
takes value k with probability proportional to S{Fn n \k) (that is, X^ is 
the number of classes in a uniformly chosen partition of F% into non-empty 
independent sets), then Xn^ is asymptotically normal, meaning that suitably 
normalized it tends in distribution to the standard normal. This generalizes 
a seminal result of Harper on the ordinary Stirling numbers. Along the way 
we give recurrences for calculating the generating functions of the sequences 
(S(F£, k))k>o, show that these functions have all real zeroes, and exhibit three 
different interlacing patterns between the zeroes of pairs of consecutive gener- 
ating functions. 

We next consider Stirling numbers of cycles. We establish asymptotic 
normality for the number of classes in a uniformly chosen partition of C n (the 
cycle on n vertices) into non-empty independent sets. We give a recurrence 
for calculating the generating function of the sequence (S(C n , k))k>o, and use 
this to give a direct proof of a log-concavity result that had previously only 
been arrived at in a very indirect way. 



* dgalvinl@nd.edu; Department of Mathematics, University of Notre Dame, Notre Dame IN 
46556. Supported in part by National Security Agency grant H98230-10-1-0364. 

^ dtrongth@nd.edu; University of Notre Dame, Notre Dame IN 46556. Supported by Glynn 
Family Honors Program. 



1 



1 Introduction 



Let G = (V, E) be a (simple, finite, loopless) graph on n vertices. An independent 
set in G is a subset of the vertices, no two of which are adjacent. For each integer 
k we set 

S(G, k) = ({partitions of V into k non-empty independent sets}| . 

We may interpret S(G, k) as the number of proper A;-colorings of G that use all 
k colors, with two colorings identified if they are identical up to the names of the 
colors. As far as we can discover, this parameter was first explicitly considered by 
Tomescu [TS]. When G = E n , the graph with no edges, S(G,k) is just the Stirling 
number of the second kind S(n, k), the number of partitions of a set of size n into 
exactly k non-empty classes. It is for this reason that in [8j, Duncan and Peele refer 
to S(G, k) as a graphical Stirling number. Motivated by the connection to colorings, 
Goldman, Joichi and White [11] refer to the sequence (S(G, k))k>o as the chromatic 
vector of G. 

Recall that associated with each G there is a polynomial Xg( x ) of degree n, the 
chromatic polynomial of G, whose value at each positive integer x is the number 
of proper x-colorings of G, that is, the number of functions / : V — > {1, . . . ,x} 
satisfying f(u) ^ f(v) whenever uv G E. The chromatic polynomial of G determines 
the sequence {S(G, k)) k > , and vice- versa. On the one hand inclusion-exclusion gives 

s ( fi! .*) = TfE(- l ) i (*)xo(*-o, (i) 

while on the other hand 

Xg(x) = ^2S(G,k)x (k) 

where x^) = x(x — 1) . . . (x — k + 1). Indeed, given a palette of x colors, for each 
k there are S(G, k) ways to partition the vertex set into k non-empty color classes, 
and ways to assign colors the classes. 

Note that S(G, k) = whenever k < x{G) (where x{G) is the chromatic number 
of G, the least positive integer x for which xg(%) > 0) and also whenever k > n. 
Also, S(G, k) is non-zero at k — x{G) an d takes value 1 at k = n. It follows that 
the generating function of the sequence (S(G, k))k>o, 

a(G,x) = J2S{G, k)x k 

k>0 

is a monic polynomial of degree n, and has a zero with multiplicity x(G) at 0. We 
refer to o~(G, x) as the Stirling polynomial of G. 
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The Stirling polynomial (arrived at via an alternate definition) was first con- 
sidered by Korfhage [13], who established some basic properties and calculated it 
explicitly for some particular families of graphs (including forests). It was then ex- 
tensively studied by Brenti [2] and Brenti, Royle and Wagner jl], who established 
large classes of graphs (including cycles) for which the Stirling polynomial has all 
real zeroes. Graphical Stirling numbers (in various guises) for some particular fam- 
ilies of graphs have also been considered by Farrell and Whitehead [9] , Duncan and 
Peele [8], Mohr and Porter [T5], Goldman, Joichi and White [UJ, and Munagi [16J. 

There are some important consequences of the generating function of a sequence 
(ak)k=o °f non-negative terms having all real zeroes. A theorem of Newton (see for 
example [T7J page 504]) implies that such a sequence satisfies the ultra log-concavity 
relation 



for k = 1, ... ,n — 1, which in turn implies the usual log-concavity relation a\ > 
Ofc-iOfc+i. If the sequence has no internal zeros (i.e., no k with au-i ^ 0, = and 
dk+i 7^ 0) then log-concavity in turn implies that the sequence is unimodal, that is, 
that there is < k < n with do < . . . < a& > . . . > a n . 

Another consequence of real zeroes is that if A is a random variable supported 
on {0, . . . , n} that takes value k with probability proportional to a^, then A may 
be realized as the sum of n independent Bernoulli random variables, leading to the 
possibility of a central limit theorem and asymptotic normality. 

Definition 1.1. A sequence (A n ) n > c is said to be asymptotically normal if, for each 
x G R, we have 



with the convergence uniform in x. 

In other words, a sequence is asymptotically normal if, suitably normalized, it 
tends in distribution to the standard normal. Note that if counts the number of 
members of size k of some family, then A measures the size of a uniformly selected 
member of the family. 

A seminal result in this area is due to Harper |12j . He proved that for each 
n, the generating function of the sequence (S(n, k))k> of the Stirling numbers of 
the second kind has all real zeroes, and used this to prove that the sequence of 
random variables (A^ tirlmg ) n > , where X^ tirling is the number of classes in a uniformly 
chosen partition of a set of size n into non-empty classes, is asymptotically normal. 
The first aim of this note is to obtain an analog of Harper's result (and, as it 
turns out, a generalization) for the Stirling numbers of forests. In what follows 
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B n = Xlfe=o S{ n > is the nth Bell number, and W(n) = logn — log log n + 0(1) is 
the Lambert W function, the unique positive solution to W(n)e w ^ = n. 

Theorem 1.2. Let (Fn ) n >o be a sequence of forests, with F% having n vertices 
and c(n) components. Let Xn be the number of classes in a uniformly chosen 
partition of the vertex set of F%^ into non-empty independent sets. There is a 
constant C > such that if c(n) < C yfnj log n for all sufficiently large n, then the 

sequence (X^ n ')„> is asymptotically normal, with 

E \ X n ') = ~ / C (n)-n „ ( 2 ) 

" +<>(*), (3) 



W(n) \logn 



and 
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\T»T(Y< n )\ - ^»>0 i ) B n+l-i ( Ei>0 ( ( j )B n -i . 

( " ) ~ L,f l t')8.-,lL,("r , )s.-,- 1 ( ) 

+o( c m. ( 5) 



W(n)(W(n) + l) \\ogn 

Conjecture 1.3. The sequence (Xn^) n >o is asymptotically normal for all 1 < 
c(n) < n. 

The denominator of (jl]) turns out to be ^ fc>0 S(Fn , k), the total number of 
partitions of F% into non-empty independent sets. For general G, the quantity 
J2k>o ^(Cfj k) is referred to in [5] as the 5eZZ number of G. For G = £' n , the Bell 
number is just the ordinary Bell number B n . 

We begin the proof of Theorem ll.2l by deriving an explicit expression for S(F£, k) 
in terms of ordinary Stirling numbers. The same expression was obtain by Korfhage 
[13] using induction; our direct derivation is based on the fact that the chromatic 
polynomial of a forest depends only on how many vertices and components it has. 
This expression allows us to obtain a number of recurrences for a(F£,x), which 
allow us to establish that o~(F£, x) has all real zeros for all n and c, (a result already 
implicit in [2J) and that moreover the zeroes of the a(i^,x)'s exhibit a number of 
nice interlacing patterns, none of which have been previously observed. 

Definition 1.4. Given two reals polynomials f and g, we say that f -< g if all of 

the following conditions hold: 

1. the zeroes of both f and g are all real and non-positive, and all negative zeroes 
of both have multiplicity one; 
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2. the number of negative zeroes of g is either the same or one greater than the 
number of negative zeroes of f ; and 



3. if Xi, . . . , x m are the negative zeroes of g and y 1; . . . , y n are the negative zeroes 
of f , both written in decreasing order, then 

xi > yi > x 2 > y 2 > . . . . 

Note that this is not the same as saying that the roots of / and g interlace in 
the normal sense; we introduce the notation / -< g to deal with the fact that most 
of the polynomials we will encounter have multiple zeroes at 0. 

Theorem 1.5. Fix c > 1 and n > c. We have the following interlacing relations 
between the zeroes of the c(F^, x) 's. 

1. a(F^ +1 ,x) -< a(F^x) 

2. cr(F°, x) -< a{F£ +l , x) for alln>c+l 

3. *(F%,x)**(F£\,x) 
I o-(F° +1 ,x)^a(F£lx) 

5. cr(F^Xl, x) -< <r{F% +1 , x) for alln>c+l 

a(F?,x) 

I \ 
a(F^x) <- <t(F%,x) 

t \ I \ 
a(Flx) -> a(Fi,x) <- a(Fi,x) 

t \ t \ I \ 
a(Flx) -> a(Flx) -> a(Flx) <- a(F 4 4 ,x) 

The interlacing patterns between the cr(F^, x)'s. 
An arrow from P to Q (with P at the head) indicates that P -< Q. 

We prove Theorem 11.51 in Section [2j where we also derive (J2]) and (jl]). Known 
asymptotic estimates for Bell numbers, together with a general theory of asymptotic 
normality (which we briefly discuss in Section HJ), then allows us to complete the 
proof of Theorem 11.21 The details are in Section [5j 

The second aim of this note is to obtain an analog of Harper's result for the 
Stirling numbers of cycles. 
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Theorem 1.6. Let X^ ycle be the number of classes in a uniformly chosen partition 
of the vertex set of the cycle C n on n vertices into non-empty independent sets. The 
sequence (X^ yclc ) n > is asymptotically normal, with 

E{XT le ) = ^-i'^' 1 (6) 



n 



Oi- (7) 



and 



W{n) \logn 



Var(X n ) = ==— ^ ~, - 1 (8) 



E J>0 (-l)^n-i-, lEi>o(- 1 ) i A 



n-l- 



+ °(^)- (9) 



W(n)(W(n) + l) \\ogn, 

In the course of the proof of Theorem II .21 we will see that the Bell number of the 
cycle C n is J2i>o(~ l) l B n -i-i- This expression also occurs in [7], where it shown to 
be the number of partitions of a set of size n into non-empty classes each of size at 
least 2. One of the main concerns of [7] is the quantity S*(n, k), the number of such 
partitions with exactly k classes, and it is shown that the random variable taking 
value k with probability proportional to S*(n, k) is asymptotically normal. We do 
not see any connection, though, between this question and the question of graphical 
Stirling numbers of cycles, other than the coincidence between the Bell number of 
C n and the total number of singleton-free partitions of [n]. 

As with Theorem 11.21 we begin by deriving a recurrence for a(C n ,x). Ideally 
we would like to use this recurrence to give a direct proof of the fact that for each 
n > 3, cr(C„,x) has all real zeroes. This was already established in jl] (after having 
been conjectured by Brenti in [2]). The methods of jl] are circuitous, involving 
results of Wagner on partition polynomials [19], a detailed study of an operator that 
converts the Stirling polynomial into the chromatic polynomial, and facts about 
the locations in the complex plane of the zeroes of the chromatic polynomial of C n . 
Unfortunately, we cannot see a direct proof. We can, however, use our recurrence 
for a(C n , x) to give a very direct proof of a strong consequence of real zeroes, namely 
ultra log-concavity of the sequence of coefficients. 

Theorem 1.7. For each n > 3 and k satisfying 1 < k < n — 1, we have 

\ 2 / 

S{C n , k) \ ( S{C n , k 




J ~ V (A) 

Moreover for k > xiPn) the inequalities above are strict. 

We prove this (together with ([6]) and in Section [3j while the remainder of 
the proof of Theorem 11.61 is in Section [5] 
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2 Forests 



Let F° be a forest on n vertices with c components (throughout this section we 
will assume that c > 1 and that n > c). It is known that a(F^,x) has all real 
zeroes; this follows from a general result of Brenti [2j Theorem 3.20], together with 
Harper's result [12] that cr(E n ,x) has all real roots. In this section we derive an 
explicit expression for S(F£,k), which allows us to obtain a number of recurrence 
relations for cr(F^, x). These recurrences allow us to give a direct proof that cr(F^, x) 
has all real zeroes, and moreover allows us to observe patterns between the zeroes of 
the Stirling polynomials of different forests (Theorem 11.51) . Our explicit expression 
for S(F£, k) also allows us to compute the mean and variance of X° (the random 
variable defined in Theorem II. 2p , a key step in the proof of Theorem 11.21 

The chromatic polynomial of F£ depends only on n and c, and not on F£ itself 
(specifically, if the components of F£ have ai, . . . , a c vertices then 

c 

Xf^x) = J] {x(x - I)*" 1 ) = x\x - l)- c ), 

and so S(F°, k) also depends only on n and c. This allows us to chose a convenient 
forest to facilitate the calculation of S(F£, k). Let F£ consist of a star onn-c-1 
vertices together with c — 1 isolated vertices. We have 

S(F n c , k) = (° T X ) S(n - 1 - 2, k - 1). (10) 

Indeed, in a partition counted by S(F°,h), the vertex at the center of the star 
cannot be in the same class as any of the leaves of the star, but there are no other 
restrictions. We get a valid partition of the vertex set into k classes by first choosing 
an arbitrary subset of size i (i G {0, . . . , c — 1}) of the isolated vertices to be in the 
same class as the center of the star, and then choosing an arbitrary partition of the 
remaining n — l — i vertices into k — 1 classes. Note that for c = 1 (the case of trees), 
the right-hand side of ffTUj) reduces to 5" (n — 1, k — 1), a fact observed in [Sj. 
By ffTUl) we have 

a(F:,x) = ^W^W-l-a-l) 

k>0 i>0 ^ ' 




where S n (x) = 'Y^ k S{n,k)x k is the generating function of the ordinary Stirling 
numbers of the second kind. This relation was first observed by Korfhage [13J, 
where an inductive proof was given. 
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It is well known that the polynomials S n (x) satisfy the recurrence 



S n ( x ) — ( x + %D)S n -i(x) 



(11) 



for n > 2 (where D is differentiation with respect to x), with the initial condition 
Si(x) = x. It follows by linearity that the polynomials a(F£,x) satisfy 



for all n > c. The initial condition for this recurrence is obtained by considering 
any forest F c c on c vertices with c components. There is only one such, the empty 
graph E c on c vertices, and so a(F^,x) = S c (x). 

We are now in a position to establish that cr(F£, x) has all real zeroes, and 
moreover that <r(F c c +1 ,x) -< <r(F c c ,x) and <r(F£,x) -< a(F£ +1 ,x) for all n > c + 1 
(items 1 and 2 of Theorem ll.5p . We begin with the first of these statements. 

From [12] we know that a(F^,x)/x is a polynomial of degree c — 1 which is 
positive at and whose roots r 1; . . . , r c _i are all real and satisfy > r x > . . . > r c _i. 
Since x(F^ +1 ,x) = 2, we have that <r(F c c +1 ,x) has a zero of multiplicity 2 at 0. We 
now consider the sign of cr(F c c +1 , x) at x — T{ for each i. By fTl2|) we have 



The first term on the right-hand side above is 0, and the second has the same sign 
as that of the derivative of cr(F c c , x)/x at x = fj. If follows that for odd i, <r(F c c +1 , r*) 
is positive, while for even % it is negative. We conclude that <r(F c c +1 ,x) has zeroes 
between r i+1 and T{ for each i G {1, . . . , c — 2}. If c is odd then lim^^oo <r(F c c +1 , x) = 
oo and so, since <r(F c c +1 , r c _x) is negative, there is a zero of <r(F c c +1 , x) below r c _i, and 
the same conclusion can be reached by similar reasoning if c is even. This accounts 
for c + 1 zeroes of <t(F c c +1 , x), and since it is a polynomial of degree c + 1 there are 
no more. It follows that <r(F c c +1 ,a;) -< o"(F c c , x). 

For item 2 of Theorem 11.51 we prove two statements by a parallel induction on n. 
First, for each n > c + 1, the polynomial a(F£,x) has all real non-positive zeroes, 
has a zero of multiplicity 2 at 0, and otherwise has distinct zeroes. Second, for 
n > c+ 1, we have cr(F^,x) -< cr(F£ +1 ,x). The base case (n = c + 1) for the first 
statement has been established in the previous paragraph. For each n > c + 1 we 
simultaneously establish the second statement for n and the first for n + 1 using an 
argument very similar to that presented in the previous paragraph, the details of 
which we leave to the reader. 

The recurrence (|12p relates the Stirling polynomials of forests with differing 
numbers of vertices on the same number of components. There is also a recurrence 




(12) 
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relating the Stirling polynomials forests on the same number of vertices with different 
numbers of components, namely 



a(F:Xlx) =a(F: +1 ,x) + a(F: 7 x) (13) 

To see this, note that 



s(K, *) = E( C i 1 V (n " 1 " i '*" 1) = ^(<- 



S(n — i,k — 1). 



i>0 v ' i>0 

and so, using Pascal's identity, 

s(F: +1 ,k)+s(F:,k) = e((-~ 1 1 ) + ( c 7 1 )) 

(n — i, k — 1) 



S(n — i, k — 1) 



i>0 v 7 



Using (Tl3|) together with items 1 and 2 of Theorem 11.51 we can easily establish 
items 3, 4 and 5 of the theorem. We begin by considering (fl3|) when n = c. If 
ri > r 2 . . . > r c _x are the negative zeroes of <r(F c c ,x) and si > s 2 > . . . > s c _i the 
negative zeroes of cr(F c c +1 , x), we know that 

> ri > si > . . . > r c _! > s c _i, 

and from (fl3|) we get that a(F^{,x) is positive at x = ri and x = Sj for all odd 
i, and negative at x = rj and x = for all even i. Since a(F^f,x) has a zero of 
multiplicity 1 at and has positive derivative at 0, it follows that a(F^l,x) must 
have a zero between n and 0, as well as between r i+ i and Sj for i e {1, . . . , c — 2}. By 
considering lim^-oo u{F^ 1 , x) we also find that there is a zero below s c _i. This 
accounts for all c + 1 zeroes, and we conclude that o"(F c c +1 ,x) -< a(F^,x) and 
cr(F£, x) -< a(F^i, x). The verification of all the remaining statements in Theorem 
11.51 is very similar, and we leave the details to the reader. 

We now derive expressions for the mean and variance of the number of 
classes in a uniformly chosen partition of F£ into non-empty independent sets, in 
terms of Bell numbers. We begin by noting that 

V, >n kS{F< k) 

£(A '- ) = few (14) 

Using (TIP]) we see that the denominator of (TT4"I) is X]«>o (^T 1 )-^™- 1 -* (this * s ^ e 
Bell number of F^, the total number of partitions into non-empty independent sets). 
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Using the recurrence S(n + l,k) = kS(n, k) + S(n, k — 1) (valid for (n, k) ^ (0, 0), 
with initial conditions S(0, 0) = 1 and S(n, 0) = S(0, k) = for n, k > 0), we get 

kS{n - 1 — i,k — 1) = S(n —i,k — l) — S(n — 1 — i,k — 2) + S(n - 1 — i,k — 1), 

and so, again using fflUj) . the numerator of f lT4|) is ^i>o (^I 1 )^™-*- This gives (T5]). 
The basic recurrence for the Stirling numbers also gives 

k 2 S{n - 1 - i, k - 1) = S(n + l-i,k-l)-2S(n-i,k-2) + 2S(n-i,k-l) 

+S(n - 1 - i, k - 3) - 3S(n - 1 - i, A; - 2) + S{n - 1 - i, k - 1). 

From this it follows that 



and so 



E((X C ) 2 ) = ^ * ^ (-^n+l-i P>n-\-i 



Sj>0 ( i )-^n-l-i 

from which (J3J) easily follows. 

3 Cycles 

We begin this section by deriving a recurrence for S(C n , k) that will allow us to com- 
pute the mean and variance of the number of colors in a randomly chosen coloring 
of the cycle, and which will also be the key to Theorem 11.71 

Let P n be the path on n vertices. The graphical Stirling numbers S(C n , k) satisfy 
the recurrence 

S(C n , k) = S(n - 1, k - 1) - S(C„-i, k) = S(P n , k) - 5(C n _i, k) (15) 

valid for all n > 4 and all k, with initial conditions S(Cz, 3) = 1 and S(C^, k) = for 
k 7^ 3. To see this recurrence, note that as a special case of ffl~2]) we have S(n — l,k — 
1) = 5(P n , k). Among the partitions of P n into k non-empty independent sets, those 
in which the first and last vertices of P n fall into different classes are in one-to-one 
correspondence with partitions of the vertex set of C n into k non-empty independent 
sets (so there are S(C n ,k) such). The remaining partitions (in which the first and 
last vertices of P n fall into the same class) are in one-to-one correspondence with 
partitions of the vertex set of C n _i into k non-empty independent sets (so there are 
S(C n -i, k) such). 

From (Tl5|) we get the polynomial recurrences 

a(C n ,x) = xS n „i(x) - cr(C n _i,x) (16) 
= a(P n ,x) - a(C n -!,x) 



10 



for n > 4, with initial condition a(Cs,x) = x 3 . Iterating (jl~6j) we find that 

n-1 

a(C n ,x) = x^(-l) i+1 ^n-i(^) 

i=l 
fc>0 i>0 

a formula that is valid for all n > 2 if we interpret C 2 to be a single edge. Recalling 
(TTTT) . we also have the recurrence 

a(C n , x)=x {{x + xD) { ^ {Cn ' x UX) ^j ) + (-l)V (17) 

for n > 3, with initial condition a(C 2 , x) = x 2 . 

From here we see no way of obtaining a direct proof that cr(C n , x) has all real 
zeroes, the problem being that although the operator x(x + xD)(l/x), when applied 
to a polynomial, preserves the property of having real zeroes, we cannot say the 
same with the addition of the (— l) n x 2 term. However, we do easily obtain a slightly 
weaker statement, the log-concavity of the sequence (S(C n , k)/{vj)^ =0 (Theorem 
II. 7p . Our strategy is to show that the x(x + xD){\/x) operator preserves ultra log- 
concavity, and then to deal with the (— l) n x 2 term by observing that it only impacts 
a bounded number of the ultra log-concavity relations, which can be dealt with by 
hand using our knowledge of the chromatic polynomial of the cycle. 

Let P(x) = J2l=2 a k xk satisfy n > 3 and > for each 3 < k < n. We begin by 
observing that if P(x) is ultra log-concave then so too is Q(x) = x(x+xD)(l/x)P(x), 
and that moreover for each k > 3 (or k > 2 if a 2 > 0) the ultra log-concavity 
inequalities for Q(x) are all strict. To see this, note that we have 

Q( x ) = a 2 x 2 + ^f2(ak-i + {k- l)a k )x k j + a n x n+1 . 

To show ultra log-concavity of Q(x) with strict inequalities for k > 3 (or k > 2 in 
the case a 2 > 0), we need to establish the three relations 

q 2 + 2a 3 \ ( a 2 \ ( ^3 + 3a 4 ^ 

m ) vet 1 ) A /' 

a n _i + (n - l)a„ \ ( a n \ ( a n _ 2 + (n - 2)a n _i . , , 

m J Ac::) A K) 

and, for /c = 4, . . . , n — 1, 

a k -i + {k- l)a k \ ( Qfc-2 + {k- 2)a fc _A / a fc + ka k+1 . 

m J A Jv g:d 1 
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After some algebra, ([TBI) is seen to be implied by (thought not equivalent to) a|/(n — 
1) > a 2 a4 : /(n — 2), which for n > 3 easily follows from (03/(3)) — ( a2 /0) (^/O); 
a consequence of the ultra log-concavity of P(x). Similarly ( |T9l) is implied by 
a l-i/{ n + 1 ) > 2a n _ 2 a n /n, a consequence of (a„_i/ ( n "J) > (a„_ 2 / ( n " 2 )) (a n / (")). 
This leaves ( 120]) . which certainly holds if each of 

(21) 
(22) 




and 






U+l/ 

(23) 

hold. From (a fc -i/( fe " 1 )) > (a*^/^)) OVO) we § et (Ell) (after some algebra 
it reduces to k 2 > k 2 - 1), and from (a k /(l)) 2 > (a fe _i/ '(^J) (afc+i/( fe +i)) we get 
(J22D (it reduces to (jfc - l) 2 (n - fc + l) 2 > ((Jfe - l) 2 - l)((n - fc + l) 2 - 1)). Finally 
(1231) follows (after some algebra) for all k in the given range (and more generally 
for n > k — 2) from another consequence of the ultra log-concavity of P(x), namely 
(WG^J) (a,/ 0) > (a fc _ 2 /(,! 2 )) (aWCjJ). 

We now prove Theorem 11.71 by induction on n. The base case n = 3 is trivial. If 
n > 3 is odd, then (by the inductive hypothesis and the observation of the previous 
paragraphs), x(x + xD)(a(C n -i, x)/x) is ultra log-concave with strict inequalities 
for k > 2, and, since C n _i is connected and bipartite, the coefficient of x 2 is 1. It 
easily follows that a(C n ,x) = x(x + xD)(a(C n -i,x)/x) — x 2 is ultra log-concave 
with strict inequalities for k > 3. For even n > 3, x(x + xD)(a(C n -i,x)/x) is 
also ultra log-concave with strict inequalities for k > 3, but, since C n _i is not 
bipartite, the coefficient of x 2 is 0. To conclude ultra log-concavity of a(C n ,x) = 
x(x+xD)(cr(C n -i,x)/x) +x 2 with strict inequalities for k > 3 we also need to show 
that 



2 



S(C n ,3)\ S(C n ,4) 



(3) J ' 00 

Using Xc n ( x ) = (— l) n (^ — 1) + (x — l) n and (JTJ) we get 

Xc n (3)-3 Xc „(2) 2" -(-If -3 



and 

S(C n ,i) 



S(C n ,3) 3 , 6 

Xc„(4) - 4 X c„(3) + 6 XCn (2) _ 3" - 4 ■ 2" + (-1)" + 6 
4! " 24 
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from which (EH) follows for n > 4. 



Our recurrence for S(C n , k) also allows us to compute the expectation and vari- 
ance of X^ cle , the number of classes in a uniformly chosen partition of the vertex 
set of C n into non-empty independent sets. Iterating ( |T5l) we get 

S(C n , k) = Y^(-l) l S(n - 1 - i, k - 1) 



and so 



and 



£(xr e ) 



i>0 



E,>o(-i) J E fc >o^- i-<,*-i 



cycle ^ _ Ei>o(-l)'E fc >ofc 2 g(n-l-»,fc-i: 



y y £feo(-i)'£*>o<%-i-i.*-i) ■ 

The calculations made in Section [2] to verify (j2J) and (TjJ also easily yield flSJ) and (JSJ) 



4 Asymptotic normality 



The connection between real zeros and asymptotic normality was first made by Levy 
[H], and rediscovered by Harper [12], who was considering the specific problem of 
the Stirling numbers of the second kind. It has since been presented in a more 
general form by numerous authors, including Bender [3J, Canfield [5] and Godsil 
[TU] . The following general result is from [3J. 

Theorem 4.1. Let (C n ) n >o be a sequence of sets, and let each C n come with a 
partition C n = U£ =0 C^. Set \C n \ = c n and |C*| = c k n . Suppose that for each n the 
polynomial P n (x) = Efc=o c n xk has all real zeroes. Let X n be the upper index of a 
uniformly chosen element oj of C n (that is, X n (u) = k if and only if u G C k n ). If 
Var(X„) — > oo as n — > oo then (X n ) n > is asymptotically normal. Moreover, we 
have the local limit theorem 



lim sup 



y/Vai(X n ) Pr [X n = E(X n ) + xy/Var(X, 



-x 2 /2 



27T 



0. 



The advantage of a local limit theorem is that it provides quantitative informa- 
tion about the c^'s that asymptotical normality does not; see for example [3J, [5] or 
[7] for details. 

For completeness we provide a short proof of the asymptotic normality part of 
Theorem 14.11 that gives explicit information about the rate of convergence. The 
probability generating function of X n is $ n ( x ) — Pn{x) / P n {l) . Since $ n (x) has all 
real zeroes and all non-negative coefficients, we may factor it as 



= n 



A; 



+ 



X 



1 + Xi 1 + A 
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with each Aj non-negative. It follows that we may write 

n 

where the I„ s are independent, and each X$ takes value with probability \/ (1 + 
Aj) and value 1 with probability 1/(1 + Aj), and so 

X n -E{X n ) = Y,Yi 

i=l 

where the Y^s are a collection of independent, mean random variables, with specif- 
ically Yi taking value — 1/(1 + Aj) with probability Aj/ (1 + Aj) and value Aj/(1 + Aj) 
with probability 1/(1 + Aj). Since Yj only takes values between —1 and 1, we have 
£(|K t | 3 ) < E(Y?). 

The Berry-Esseen theorem [1] now says that there is an absolute constant C > 
such that for all x G M we have 

< C 
~ y/Vax(X n y 

from which asymptotic normality follows. 



Pr [ ^=i^ < x ) _ 1 f e - x2 ' 2 dx 



5 Proofs of Theorems 1.2 and 1.6 



We begin by establishing ([3]) and (JH]). We need two preliminary estimates involving 
Bell numbers. The first is due to Canfield and Harper [6]. We have 



W{n) 



B r , 



11 



O 



W{n) 



O 



n 



logn 



n- 



(25) 



(where recall W(n) = logn — log logn + 0(1) is the Lambert W function, the unique 
positive solution to W(n)e w ^ = n). The second is due to Harper [12]. We have 



B 



n+2 



B,, 



B n +i 

~b7 



n 



W(n)(W(n) + r 



+ 



log 2 n 



(26) 



(Harper gives an error term of o(l); see for example [7] for the more refined error.) 
We need a slight strengthening of ( |25l) . For k = o(n), fl25l) gives 



S n _ fc _! W(n - k) 



B 



+ 



n—k 



n 



log(n — k) 
{n - k) 2 



Win) ^ ( A; logn 
— — + 1 



n 



n- 



(27) 
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with the same bound for B n+ k-i/ B n+ k- Here we have used the mean value theorem 
to estimate W (n) - W (n - k) , noting that W (x) = W (x) / (x ( 1 + W (x) ) ) . A similar 
argument gives 

B n+i _ ( B n \ 2 = n f 1 \ 

B n - X \B n „J W(n)(W(n) + l) + \\og 2 n) ' 1 J 

From ff27|) we get that for m G {n. — 1, n, n + 1} we have 



n I \ \ n 2 



(the second equality valid for c = o(n/v / logn)). 

It now immediately follows that E(X^) = {n/W{n)){l + 0(l/n)) (recall that we 
are in the range c = 0{^Jn/ logn)). It also follows that 

l'l / 1 \ ^ 2 



+ 



cr 



o' 1 



W(n)(W(n) + l) \\ogn 

(here using ff271) and fl28]) ). As long as c < C\Jnj log n for some sufficiently small 
C > 0, we have Var(A ra ) = w(l), and so we complete the proof of Theorem 11.21 by 
appealing to Theorem 14.11 

We now turn to Theorem 11.61 For each k G { — 1, 0, 1} we use (I27p to obtain 



Bn+k-i — B n+k ( 1 
i>0 ^ 



Bn+k-i , B n+ k_ 2 B n +k-i 



B n +k B n+ k-iB n - 



-k 



B n+ k-3B n +k-2B n +k-i , q ( nB n+k _ 4 



B n +k-2B n +k-iB n+k \ B n+ k 



n n 2 n 3 / V \ n 2 



Using and (jEJ) this gives 



E(X ^) = #s- (i + o = + o 



S„_x V \ n2 J J W(n) \logn 
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n ( 1 

W(n)(W(n) + 1) + V.logn 

the last equality from fl28|) . Since Var(X£ ycle ) = w(l), we complete the proof of 
Theorem 11.61 by appealing to Theorem 14.11 



References 

[1] A. Berry, The Accuracy of the Gaussian Approximation to the Sum of Indepen- 
dent Variates, Transactions of the American Mathematical Society 49 (1941), 
122-136. 

[2] F. Brenti, Expansions of chromatic polynomials and log-concavity, Trans. Amer. 
Math. Soc. 332 (1992), 729-756. 

[3] E. Bender, Central and local limit theorems applied to asymptotic enumeration, 
J. Comb. Theory Ser. A 15 (1973), 91-111. 

[4] F. Brenti, G. Royle and D. Wagner, Location of zeros of chromatic and related 
polynomials of graphs, Can. J. Math. 46 (1994), 55-80. 

[5] E. Canfield, Central and local limit theorems for the coefficients of polynomials 
of Binomial type, J. Comb. Theory Ser. A 23 (1977), 275-290. 

[6] E. Canfield and L. Harper, A simplified guide to large antichains in the partition 
lattice, Congr. Numer. 100 (1994), 81-88. 

[7] E. Czabarka, P. Erdos, V. Johnson, A. Kupczok and L. Szekely, Asymptotically 
normal distribution of some tree families relevant for phylogenetics, and of 
partitions without singletons, Moscow J. Comb. Numb. Theory 1 (2011). 

[8] B. Duncan and R. Peele, Bell and Stirling Numbers for Graphs, Journal of 
Integer Sequences 12 (2009), article 09.7.1. 

[9] E. Farrell and E. Whitehead, Connections between the Matching and Chromatic 
Polynomials, Internal J. Math. & Math. Sci. 15 (1992), 757-766. 

[10] C. Godsil, Matching behaviour is asymptotically normal, Combinatorica 1 
(1981), 369-376. 



16 



[11] J. Goldman, J. Joichi and D. White, Rook Theory III. Rook polynomials and 
the Chromatic structure of graphs, J. Combin. Th. Ser. B 25 (1978), 135-142. 

[12] L. Harper, Stirling Behavior is Asymptotically Normal, Ann. Math. Statist. 38 
(1967), 410-414. 

[13] R. Korfhage, cr-polynomials and graph coloring, J. Comb. Theory Ser. B 24 
(1978), 137-153. 

[14] P. Levy, Sur une propriete de la loi de Poisson relative aux petites probabilites, 
Soc. Math, de France, Comptes rendus des sceances de Vannee 1936, (1937), 
29. 

[15] A. Mohr and T. Porter, Applications of Chromatic Polynomials Involving Stir- 
ling Numbers, Journal of Combinatorial Mathematics and Combinatorial Com- 
puting 70 (2009), 57-64. 

[16] A. Munagi, /c-Complementing subsets of nonnegative integers, Internat. Jour. 
Math. 21 (2005), 215-224. 

[17] R. Stanley, Log-concave and unimodal sequences in algebra, combinatorics and 
geometry, Annals of the New York Academy of Sciences 576 (1989), 500-534. 

[18] I. Tomescu, Methods combinatoires dans le theorie des automates finis, Ph. D. 
thesis, Bucarest, 1971. 

[19] D. Wagner, The partition polynomial of a finite set system, J. Comb. Theory 
Ser. A 56 (1991), 138-159. 



17 



